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we also discuss some issues related to the apparent and event horizons of the black 
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1. Introduction 

In 1951, P. C. Vaidya presented for the first time a metric to describe the spacetime 
outside a spherically symmetric star, taking into account the radiation flux emitted by 
the star [TJ. In his work, the mass of the star is no longer considered to be constant and 
the metric is not static. The situation is described as a spherical mass surrounded by a 
finite and nonstatic envelope of radiation with radial symmetry. 

This metric has been the usual starting point for the study of the quasinormal 
modes of time dependent black holes [21 EJ EJ |6] . The charged version of the metric is 
also an usual starting point for the study of many aspects of charged black hole physics 

ei ei ina mi H2i . 

The Vaidya metric has also been used to describe spherically symmetric collapse 
and the formation of naked singularities [HI fT5] . It was also applied to the study of 
Hawking radiation and the process of black hole evaporation [161 EH (HI ESL m the 
stochastic gravity program [20J, and in recent numerical relativity investigations [2Tj . 

In the study of quasinormal modes, it is useful to write the Vaidya metric in double 
null coordinates [H [5j [22] . However, there is no general coordinate transformation from 
the usual radiation coordinates to double null coordinates [23]. So our purpose with this 
paper is to present the most general Vaidya metric, n-dimensional, with cosmological 
constant A and with time dependent mass and electrical charge in these coordinates. 

The numerical results are obtained with a generalization of the semi-analytical 
method used in [211 EE] , that allows us to construct the structure of the spacetime from 
the behavior of the outgoing null geodesies. 

The structure of the paper is as follows. In section [2] we present the derivation of 
the most general Vaidya metric in double-null coordinates. In section |3] we present a 
discussion on the structure of the spacetime and the numerical results obtained with 
the semi-analytical method. And finally in section H] we present our final discussion of 
the results. 

2. The general Vaidya metric in double- null coordinates 

The n-dimensional Vaidya metric was first discussed in [26]. It can be easily cast in n- 
dimensional radiation coordinates (v, r, 9i, . . . , Q n -2) as done, for instance, in [27]. The 
n-dimensional charged Vaidya metric in radiation coordinates, obtained originally in [8] , 

reads 

ds2 = " f 1 " (^fU + ( B _ 2,(^3)^) ) d " 2 + 2 * d " + W 

where n > 3, c = ±1, and df2^_ 2 stands for the metric of the unit (n — 2)-dimensional 
sphere, assumed here to be spanned by the angular coordinates (6*i, 6*2, ... , 9 n -2), 

n-2 / i-1 \ 

c a =E n^-k- (2) 
»=i \j=i j 



Double null formulation of the general Vaidya metric 



3 



For the case of an ingoing radial flow, c = 1 and m(v) is a monotonically increasing mass 
function in the advanced time v, while c = — 1 corresponds to an outgoing radial flow, 
with m(v) being in this case a monotonically decreasing mass function in the retarded 
time v. The constant q corresponds to the total electric charge. In principle, one can 
also consider time dependent charges q as done, for instance, in [9]. This situation will 
of course require the presence of charged null fluids and currents, whose realistic nature 
we do not address here. 

It has been known since a long time that the radiation coordinates are defective 
at the horizon [28], implying that the Vaidya metric (pQ), with or without electric 
charge, is not geodesically complete in any dimension. The radiation coordinates 
are not enough to cover the entire Vaidya spacetime. (The radiation coordinates are 
defective at horizons where v 2 — > oo). As can be seen in [28], for a 4-dimensional 
Vaidya metric with dm/du < and without electric charge (in the context of ref. [28J, 
a radiating star), the hypersurface r = 2m(oo) at v = oo in the Vaidya metric is 
analogous to the Schwarzschild hypersurface r = 2m at Schwarzschild's time coordinate 
T = +oo in the Kruskal metric. (See [10] and [29] for further discussions about possible 
analytical extensions and properties of the horizon of the Vaidya metric in the radiation 
coordinates) . 

The cross term drdv introduces extra terms in the hyperbolic equations governing 
the evolution of physical fields on spacetimes with the metric (pD). Typically, the 
double-null coordinates are far more convenient for QNM analysis. This was the main 
motivation of the series of works based on Waugh and Lake's approach [23J, where 
the problem of casting the 4-dimensional Vaidya metric in double-null coordinates was 
originally addressed. As all previous attempts to construct a general transformation 
from radiation to double-null coordinates had failed, Waugh and Lake considered the 
problem of solving Einstein's equation with spherical symmetry directly in double-null 
coordinates. The resulting equations, however, are not analytically solvable in general. 
Waugh and Lake's work was revisited in [23], where a semi-analytical approach allowing 
for general mass functions was proposed. 

More recently this semi- analytical approach was extended to the case of an n- 
dimensional Vaidya metric with cosmological constant A [25]. This approach consists 
in a qualitative study of the null-geodesics, allowing the description of light-cones and 
revealing many features of the underlying causal structure. It can also be used for more 
quantitative analyses; indeed, it has already enhanced considerably the accuracy of the 
quasinormal modes analysis of varying mass black holes [HE], and it can also be applied 
to the study of gravitational collapse [24]. 

In this section, we extend the approach proposed in |25j and derive the double- 
null formulation for the most general Vaidya metric: n- dimensional, in the presence of 
a cosmological constant, and with varying electric charge and mass. Only the main 
results are presented. The reader can get more details on the employed semi-analytical 
approach in [25] and the references cited therein. We recall that the n-dimensional 
spherically symmetric line element in double-null coordinates ( m, v, . . . , Qn—2) is given 
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by 

ds 2 = -2/(«, v)du dv + r 2 (u, v)dtt 2 n _ 2 , (3) 

where f(u,v) and r(u,v) are non vanishing smooth functions. The energy-momentum 
tensor of a unidirectional radial null-fluid in the eikonal approximation in the presence 
of an electromagnetic field F ab is given by 

T ab = ^-hk a k b + - 1 - (F ac F b c - -g ab F cd F cd ) , (4) 

87T 47T \ 4 / 

where k a is a radial null vector and h(u, v) is a smooth function characterizing the null- 
fluid radial flow. In the double null coordinates (u, v,9x, . . . , 6 , „_ 2 ), we can choose either 
k a = (1,0,0,..., 0) (flow along the u-direction) or k a = (0, 1, 0, ... , 0) (flow along the 
v-direction). Since the u and v directions are unspecified, it is not in fact necessary to 
consider flows along both directions. We will consider here, without loss of generality, 
the case of a flow along the v -direction, as done in [23] . 
Maxwell equations are given by 

1 --daV=gF ab = 4vrJ 6 , (5a) 



F^c + F ca>b + F bCia = , (56) 
where, for the metric ([3]), 

n—3 

v / ^=/r n - 2 Il(sm^r^ 2 . (6) 

i=i 

All geometrical quantities relevant to this work are listed in the appendix of [25]. The 
equations (|5aj) have the following static spherically symmetric solution 

/ 'ur i <rn 'I 



j r n-2 ' 

with all other components of the electromagnetic tensor vanishing, where q is a constant 
which one identifies as the n-dimensional electric monopole charge. This case, of course, 
corresponds to J b = 0. In order to allow for a time dependent charge q(v), one needs to 
assume the presence of a current 

ru _ 1 gfo) /OX 

- 4?r/r n-2' W 

(with all other components vanishing) which is obtained from the continuity equation 
J a . a = 0, as done, for instance, in [9]. Such a current must naturally appear, as we will 
see, in the h(u, v) function characterizing the radial flow in the energy momentum tensor 
01]). We note here that our solution ([7]) for F ab with q(v) also consistently satisfies the 
sourceless Maxwell equations (156|) . 

Einstein equations with cosmological constant A 

Rab - T^QabR = -hg ab + 87rT a6 , (9) 

imply that, for the energy-momentum tensor (j3J), 

2n t n — 4 q 2 , , 

R= A-2 -f— , 10 

n-2 n -2r 2 ( n ~ 2 y v ; 
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where fl7]) was used. The uu and m components of Einstein equations for this case read 
simply 

f r >u 

^_I> = _& _1 (12) 

where >w and >v denote, respectively, differentiation with respect to u and v as usual. 
The uv and $.g. components are, respectively, 

^ - f ~f - (n - 2)^ = -^-/ + 2^|/^, (13) 
/ r n — 2 n — 2 r 2 ( n 2 > 

2 2A 2 q 2 

- (rr tUV + (n - 3)r u r „) + (n - 3) = -r 2 + (14) 

/ n — 2 n — 2 r^ n A > 

For n ^ 3, differentiating Eq. ( fl4l) with respect to it and then inserting Eq. ( fTTT) leads 

to 

r w - 2 r,„ 2A w-1 2 q^_ _ 

f {n-2){n-\) r n-2r n ~ 3 ™' 1 ' 

(after integration with respect to it) where m(v) is an arbitrary integration function that 
we already known from [25] that can be interpreted as the mass of the solution. The 
n = 3 case must be considered separately, in an analogous way as done for q = in [25J, 
and the most important results are presented in the Appendix. Now, differentiating Eq. 
( |T4|) with respect to v and using ( fl2l) and ( IIBT) gives 

ft= _(^^_ L. _ i ijtk). (16) 

Vn — 3/ r n u \ ' n — 2r nA ) 
Eq. (TTTT) is ready to be integrated 

f = 2Br, u , (17) 

where -B(v) is another arbitrary (but nonvanishing) integration function. From (Tl6|) and 
(TTTj) . one has 



v rt — 3/ r n_2 V ' n — 2 r n_3 / 
Finally, by using (I14j) and (fl71) . Eq. ( Tl~5l) can be written as 

r = _ B / r 1 ^ ?A 2 2_ <ft±\ 

V (n-3)r™- 3 (n-2)(n-l) (n - 2)(n - 3) r 2 (™" 3 ) / ' 1 ; 

Note that f lT3|) and ( TTTJ) reproduce (ITS"]) . Einstein equations are, therefore, equivalent 
to the equations f|T7|) . f|T8|) . and f|T9|) . generalizing the previous results of [23], [21] and 
[25]. 

As already mentioned, the physical interpretation of the arbitrary integration 
functions m(v) and B(v) are the same of the q = case. Transforming from the 
double-null coordinates back to the radiation coordinates by the coordinate change 
(u,v) — » (r(w, t>), v ), the metric (jS]) will read 

ds 2 = 4Br, v dv 2 - ABdrdv + r 2 dtt 2 n _ 2 , (20) 
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where (ITT]) was explicitly used. Comparing ((T]) and (12"01 and taking ( TIT)]) into account, it 
is clear that with the choice B = ±1/2, the function m(v) indeed represents the mass of 
the n-dimensional charged solution. The coordinate transformation leading to ( 120]) also 
ensures that the Vaidya metric in radiation coordinates (JT|) and the double null metric 
()3]) constructed in this paper are (locally) isometric. It is important to stress this fact, 
given the absence of a Birkhoff theorem for non vacuum spacetimes. 

As in the q = case, the weak energy condition applied for (rjj requires, from f ]T8]) . 

that 

bU,~Z 1 - 2 ^S)<0- (2D 



n 



If there are both mass and charge variations, then m^ v and q^ v cannot be chosen 
arbitrarily (see [9] for a discussion) and must chosen satisfy the energy condition ff2T]) . 

Taking (as mentioned above) B = ±|, if we have only mass (or charge) varying 
with time, the energy condition requires that m(v) (or q(v)) be a monotonic function 
and fixes B in the following way: 

if m v > (or q v < 0) then B = — , 
if m >v < (or q >v > 0) then B = + - , 
where we consider, without loss of generality, q{y ) > 0. 



3. The spacetime structure 

The problem of constructing a double-null formulation for the general Vaidya metric 
may be stated as follows: given the mass function m(v), the electric charge function 
q(v), the cosmological constant A, and the constant B, one needs to solve Eq. ( TT9]) . 
obtaining the function r(u,v). Then, f(u,v) and h(u,v) are calculated from (fTTj) and 
(|T8]) . The arbitrary function of u appearing in the integration of ( TT9]) must be chosen 
properly [23] . so that f(u, v) given in ( ]T7j) is a non- vanishing function. Unfortunately, as 
stressed previously by Waugh and Lake[23], such a procedure is not analytically solvable 
in general. In [25], a semi-analytical procedure is introduced to attack the problem of 
solving Eqs. ( TT7]) -( TT9]) for the q = case, generalizing in this way the results of [21] 
obtained for n = 4 and A = 0. The approach, which we will not reproduce here, allows 
us to qualitatively construct conformal diagrams, identifying horizons and singularities, 
and also to evaluate specific geometric quantities. The main idea, however, is to solve 
Eq. ( 1T9]) numerically as an initial value in v, for constant u. In other words, we obtain 
numerically r(u,v) for u constant, starting with a initial condition 

r(u,v ) = F(u), (22) 

where we must have F'(u) ^ 0, as can be seen from (|T71) . In analogy with the flat 
spacetime case, we choose here F(u) = — | . Since the lines of constant u (or v) are null 
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Figure 1. Example of u-constant null- geodesies (dotted lines) obtained from Eq. (fT9|) 
for the usual Reissner-Nordstrom case, corresponding to B = —1/2, m = 0.5, q = 0.4, 
n = 4, A = 0, and taking r(u,v — 0) — —u/2. (See [21]). Here r + and r_ (solid 
lines) are defined as usual as r± — m ± y m 2 — g 2 . In the exterior region (r > r + ), 
the constant u null geodesies reach Z + , while in the interior region they are confined, 
giving origin to the typical black hole causal structure. 

geodesies for any metric in double-null coordinates, knowing r(u,v) for u constant is 
enough, for instance, to construct the causal conformal diagrams. Figure ([[]) depicts a 
simple example, corresponding to the usual Reissner-Nordstrom solution. 

Figures |2] and [3] present the behavior of the ^-constant null-geodesics in two different 
time dependent cases: increasing charge function q(v) and increasing mass function 
m(v). Charge and mass variations produce opposite results for the horizons. Also, r_ 
and r + are now time dependent, and the event horizon r^ is no longer coincident with 
the apparent horizon r + . 

The case presented in Fig. 2 deserves a more detailed analysis. The charge increase 
is analogous to the mass evaporation studied in [2U 125] ■ 

Following our discussion in Section [21 when g „ > we must have B = +|, in order 
to satisfy the weak energy condition f|2ip . Now there is a subtlety regarding the sign 
of B. We can see from Eq. ( JT7]) that the sign of / depends on the signs of B and r u 
(which is of negative sign with our choice of F(u)). Therefore, for B = — |, we have 
/ > and d u + d v is timelike and d v — d u is spacelike. However, for B — +|, / has the 
opposite sign and the timelike and spacelike directions are now exchanged. 
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r 

Figure 2. Example of it-constant null-geodesics (dotted lines) obtained from Eq. (|19|) . 
with the same parameter values used in FigureQ]but now with a time dependent charge 
function given by 2q(v) = (qf + qi) + (qf — qi) tanh p(v — v m ), with qi = 0.4, qf = 0.48, 
p = 4.0 and v m = 1.5. We can see in this case that the r± (solid lines) are no longer 
constant, and the event horizon r/j (dashed line) no longer coincides with r+. 



The transformation (u, v) — > (v,—u) restores the temporal and spatial directions. 
Under this transformation, Eq. (|T9l) (with B — +~) becomes 

1 / 2m(-u) 2A 2 2 q 2 (-u)\ 

~ T ' U ~~2 \ ~ (n-3)r™- 3 ~ (n-2)(n- if + (n-2)(n-3) r 2 ("" 3 ) J ' 

which is formally identical to a case with B = — ~, mass function m(—u) and 
charge function q(—u). Therefore, the transformed equation results in a situation with 
decreasing charge, that is, a time reversal of the original situation with increasing charge. 
Stated in a different way, this means that the case with increasing charge and B = +| 
must be interpreted as the time reversal of the case with decreasing charge and B = — | 
(with both cases satisfying the weak energy condition). 

However, in order to describe an actual charge increase (equivalent to an 
evaporation), we need to violate the weak energy condition, since there are no classical 
processes that can lead to black hole evaporation. In order to do that, we deliberately 
choose B — — | together with q jV > 0. The weak energy condition is violated, and the 
resulting evaporation process is shown in Figure 2. 
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Figure 3. Same as Figure EJ but this time with q — 0.4 and 2m(v) — (to/ 4- rrij) + 
(to/ — TOi)tanhp(u — V m ), with = 0.5, m/ = 0.56, p — 4.0 and u m = 1.5. 

In Figure H] we show an example of a higher dimensional spacetime with n = 5. The 
results are qualitatively the same as for the 4- dimensional cases we presented in Figs. [2] 
and El In Figure [5] we have a qualitatively different behavior, due to the appearance of 
a cosmological horizon (A ^ 0). 

The calculation of the apparent horizons r ± (v) shown in Figures [TJt5] follows the 
semi-analytic method described in refs. [211 [25] . We can see that the curves defined by 
the vanishing of the r.h.s of Eq. ( fl9l) always describe the frontier between two regions 
of the (v,r) plane where the solutions of Eq. ( f!9l) have distinct behaviors. For r„ < 0, 
the null geodesies approach the singularity, whereas for r v > 0, the null geodesies tend 
to escape from the singularity. 

The determination of the event horizon is done numerically by inspection of 
the initial value r(u, Vq). The event horizon is found as the last geodesic that escapes 
towards infinity and does not fall into the singularity. Note that in the case of increasing 
charge (or, equivalently, decreasing mass) there are null geodesies that escape towards 
infinity even though they were initially inside the apparent horizon r+. 
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Figure 4. Same as Figure G2 but this time exploring a 5-dimensional solution 
(ri = 5), with q — 0.2, m, = 0.5 and nif = 0.65. The apparent horizons are now 

and therefore we must have y/3m > 2q. The 
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behavior is qualitatively the same as for the n = 4 case. 



4. Final discussion 

We have presented a formulation in double-null coordinates of the most general Vaidya 
metric: n-dimensional, with varying mass and/or charge and cosmological constant A. 

By exploring the numerical solutions of Eq. f lT9|) . we were able to highlight some 
interesting features of the behavior of time-dependent horizons in multiple-horizon 
spacetimes. The w-constant geodesies can be used to track the time dependent event 
and Cauchy horizons, that no longer coincide with r + and r_. 

The formulation presented here for the Einstein Eqs. f ll7p - (TT9|) was recently used 
in a quasinormal mode analysis of the Vaidya metric [22], and provided the framework 
needed to obtain the quasinormal frequencies with sufficient accuracy to verify their 
nonstationary behavior. 
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r 

Figure 5. Same as Figure [31 but this time allowing a non-zero cosmological constant 
A = 0.3, with q = 0.4, rrii — 0.5, rrif = 0.65, p — 1.0 and v m — 15. The striking feature 
in this case is the existence of three time-dependent apparent horizons r_,r + and r c , 
located at the positive roots of Eq. (|19l) , which is a 4th order polynomial in this case. 



Appendix 

Here we present a short discussion and the results for the Einstein equations (TTTT)-(TT4l) 
for the n = 3 case, generalizing the discussion of |25| . For n = 3, Eq. (TTTj) is still valid 
and can be integrated to give 

f = 2Br tU , (A.l) 

same as Eq. (ITT]) . 

Taking n — 3, Eq. ( THE]) reads 

rr ul ,-(Ar 2 + g 2 )/ = 0, (A.2) 

and we can use Eq. (lA.lj) to integrate Eq. flA.2j) and obtain 

r>v = -B(-m - Ar 2 - 2g 2 In r) , (A.3) 

which is the n = 3 version of Eq. (1191) . Here m(v) is an integration function that has 
the same interpretation as before, that is, the mass of the solution. Compare (1A.3j) with 
the charged BTZ black hole [30] . 



Double null formulation of the general Vaidya metric 



12 



From Eq. (1121) with n = 3, we have 

/i = ^%- — , ( A -4) 
r f r 

and using Eqs. (1A.3I) and flA.lj) to get 



/, = 25 Ar+^J/, (A.5) 

we obtain the n = 3 version of Eq. (TX2j) : 

n /m„ 41nr \ , , . 

h = -B\-± + ——qqA. (A.6) 

We also note here that the n = 3 version of Eq. ( {TBI) , 

f,uf,v f,uv r 



-2Af, (A.7) 



/ 2 / 

together with Eq. ( lA.lj) still reproduce Eq. ( [15]) with n = 3. 

So finally we can conclude that the Einstein equations (llll)- (jT4]) . which in the n > 4 
case are equivalent to Eqs. ([T7]) - (jT9]) . are equivalent in the n = 3 case to Eqs. (1A. 1 j) . 
flA~3|) and ftKM . 
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